本 章 主 要 内 容 和 学 习 目 的 


首先 介绍 计算 机 仿真 技术 和 仿真 建 模 方 法 的 基本 概念 ， 以 
便 对 建 模 和 仿真 有 个 初步 和 整体 的 认识 ;然后 对 Simulink 进 
行 简单 介绍 ， 并 以 一 个 简单 例子 进行 引导 ， 最 后 介绍 Simulink 
的 工作 原理 ， 为 后 续 的 深入 营 握 Simulink 打下 基础 。 


@ 动态 系统 的 计算 机 仿真 
@ 仿真 三 要 素 

@ Simulink 与 建 模 仿真 
@ Simulink 的 安装 

@ 创建 一 个 简单 模型 

@ 横 型 基本 结构 


1.1 动态 系统 的 计算 机 仿 其 


1.1.1 系统 与 模型 
1. 系统 

系统 只 指 具 有 某 些 特定 功能 、 相 互联 系 、 相 互 作用 的 
元 素 的 集合 。 这 里 的 系统 是 指 广义 上 的 系统 ， 泛 指 上 自然 界 
的 一 切 现象 与 过 程 ， 例 如 工程 系统 如 控制 系统 、 通 讯 系 统 
等 ， 非 工程 系统 如 股市 系统 、 交 通 系 统 、 生 物 系 统 等 。 


2. 系统 模型 


系统 模型 是 对 实际 系统 的 一 种 抽象 ， 是 对 系统 本 质 《 或 
是 系统 的 某 种 特性 ) 的 一 种 描述 。 模 型 具有 与 系统 相似 的 特 
性 。 好 的 模型 能 够 反映 实际 系统 的 主要 特征 和 运动 规律 。 


模型 可 以 分 为 实体 模型 和 数学 模型 。 


实体 模型 义 称 物理 效应 模型 ， 是 根据 系统 之 间 的 相似 性 而 建 
立 起 来 的 物理 模型 ， 如 建筑 模型 等 。 

数学 模型 包括 原始 系统 数学 模型 和 仿真 系统 数学 模型 。 原 始 
系统 数学 模型 是 对 系统 的 原始 数学 描述 。 仿 真 系统 数学 模型 
是 一 种 适合 于 在 计算 机 上 演算 的 模型 ， 主 要 是 指 根据 计算 机 
的 运 鼻 特 点 、 仿 真 方式 、 计 和 方法、 精度 要 求 将 原始 系统 数 
学 模型 转换 为 计算 机 程序 。 

和 动态 系统 模型 
连续 系统 模型 离 知 系统 模型 


代数 方程 一 分 布 多 数 | 参数 差分 方程 





1.1.2 计算 机 仿真 


1. 仿真 的 概念 

仿真 是 以 相似 性 原理 、 控 制 论 、 信 息 技 术 及 相关 和 领域 的 
有 关 知 识 为 基础 ， 以 计算 机 和 各 种 专用 物理 设备 为 工具 ， 借 
助 系统 模型 对 真实 系统 进行 试验 的 一 门 综合 性 技术 。 


2 

(1) 实物 仿真 : 又 称 物 理 仿 真 。 是 指 研制 某 些 实体 模型 ， 
使 之 能 够 重 现 原 系统 的 各 种 状态 。 早 期 的 仿真 大 多 属 
ia 

优点 直观， 形象 ， 至 今 仍 然 广泛 应 用 。 

缺点 : 投资 巨大 、 周 期 长 ， 难 于 改变 参数 ， 灵 活性 差 。 


(2) 数学 仿真 : 是 用 数学 语言 去 描述 一 个 系统 ， 并 编制 程 
序 在 计算 机 上 对 实际 系统 进行 研究 的 过 程 。 

优点 : 灵活 性 高 ， 便 于 改变 系统 结构 和 参数 ， 效 率 高 
(可 以 在 很 短 时 间 内 完成 实际 系统 很 长 时 间 的 
动态 演变 过 程 ) ， 重 复 性 好 

缺点 :对 某 些 复杂 系统 可 能 很 难 用 数学 模型 来 表达 ， 或 
者 难以 建立 其 精确 模型 ， 或 者 由 于 数学 模型 过 
于 复杂 而 难以 求解 


43) 半 实 物 仿真 : 叉 称 数学 物理 仿真 或 者 混合 仿真 。 为 了 
提高 仿真 的 可 信和 度 或 者 针对 一 些 难以 建 模 的 实体 ， 在 
系统 研究 中 往往 把 数学 模型 、 物 理 模 型 和 实体 结合 
来 组 成 一 个 复杂 的 仿真 系统 ， 这 种 在 仿真 环节 中 存在 
实体 的 仿真 称 为 半 物 理 仿 真 或 者 半 物 理 仿真 ， 如 飞机 
半 实 物 仿真 等 。 


3. 计算 机 仿真 


计算 机 仿真 是 在 研究 系统 过 程 中 根据 相似 性 原理 ， 利 用 计 
算 机 来 盘 真 模拟 研究 系统 。 研 究 对 象 可 以 是 实际 的 系统 ， 也 可 
以 是 设想 中 的 系统 。 在 没有 计算 机 以 前 ， 仿 真 都 是 利用 实物 或 
者 它 的 物理 模型 来 进行 研究 的 ， 即 物理 仿真 。 物 理 仿真 的 优 所 
是 直接 、 形 象 、 可 信 ， 缺 点 是 模型 受 限 、 易 破坏 、 难 以 重用 。 
计算 机 仿真 可 以 用 于 研制 产品 或 设计 系统 的 全 过 程 ， 包 括 方案 
论证 、 技 术 指标 确定 、 设 计 分 析 、 故 障 处 理 等 各 个 阶段 。 如 训 
练 飞行 员 、 宇 航 员 的 方针 工作 侣 和 仿真 机 舱 等 。 


1.2 仿真 的 三 要 素 


计算 机 仿真 的 三 个 基本 要 素 是 系统 、 模 型 和 计算 机 ， 联 
系 看 它们 的 三 项 基本 活动 是 模型 建立 、 仿 真 模型 建立 〈 又 称 
二 次 建 模 ) 和 仿真 试验 。 


数学 仿真 有 用 数学 模型 ， 用 数学 语言 对 系统 的 特性 进行 摘 
述 ， 其 工作 过 程 是 : 
@ 建立 系统 的 数学 模型 


@ 建立 系统 仿真 模型 ， 即 设计 算法 ， 并 转化 为 计算 机 程序 ， 使 
系统 的 数学 模型 能 为 计算 机 所 接受 并 能 在 计算 机 上 运行 

@ 运行 仿真 模型 ， 进 行 仿真 试验 ， 再 根据 仿真 试验 的 结 宋 进 
一 步 修正 系统 的 数学 模型 和 仿真 模型 。 


1.3 Simulink 与 建 模仿 真 


1.3.1 Simulink 


Simulink 是 一 种 用 来 实现 计算 机 仿真 的 软件 工具 。 它 是 
MATLAB 的 一 个 附加 组 件 ， 可 用 于 实现 各 种 动态 系统 〈 包 
括 连续 系统 、 离 散 系统 和 混合 系统 ) 的 建 模 、 分 析 和 仿真 。 


特点 :; 易学 易 用 ， 能 够 依托 MATLAB 提 供 的 丰富 的 仿真 资源 


1.3.2 Simulink 的 应 用 领域 


(1) 通讯 与 卫星 系统 ; (2) 航空 航天 系统 ; 
(3) 生物 系统 ; (4) 船舶 系统 ; 
(S$) 汽车 系统 ; (6) 金融 系统 ; 


(7) 控制 系统 。 


在 Matiab 命 令 


右 图 所 示 的 模型 用 来 
模拟 双 质 量 一 弹 壬 系 
统 在 光 请 平面 上 受 一 
个 周期 力 情况 下 的 运 
动 状 态 ， 其 中 周期 力 
只 作用 在 左边 的 质量 
块 上 。 
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此 模型 中 使 用 了 状态 判断 和 LQR 控 制 。 


模型 中 还 有 几 个 “隐藏 了 真实 身份 ”的 子 系统 ， 如 图 中 的 
Inputs&Sensors 模块 和 State estimator 模块 。 双击 后 可 看 到 
它们 的 “真实 面目 ”。 
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个 Xx = Ax+Bu @@ 
__ Double ee ， ee m 
y = Cx+Du 0 活 (Double clic ore info) 
In1 Ut To start a 人 simulation， use the 
State select "Simulatio 人 1-dow 


estimator 


运行 菜单 选项 【Simulation>Start】， 则 屏幕 上 出 现 双 质量 一 
弹 敬 系统 运动 状态 的 动画 图 形 。 


模型 中 的 Actural Position 模 块 和 Estimated _ Position 模块 用 来 
显示 在 一 个 周期 内 的 左边 质量 块 的 位 置 轨迹 。 


1.4 Simulink 的 安装 


系统 要 求 : 
奔腾 100 以 上 CPU，16MB 以 上 内 存 ，Windows 9x 或 Windows NT 


安 儿 ; 
随 MATLAB 安 闭 或 单独 安 钱 。 


有 启动 Simulink: 


首先 启动 MATLAB ， 然 后 在 MATLAB 窗口 中 单 击 上 面 的 
Simulink 按 钮 或 在 命令 窗口 中 输入 simulink。 


两 个 例子 


X(1) = Sin(1) X(t) 三 一 COS(t 和 十 人 
x(0)=0 C 为 常数 利用 初始 条 什 CC 一 ] 


X(1) 三 一 cos( 人 +1] ee 


创建 模型 步骤 1 xD)= -cos(D)+1 


X(U) = 
步骤 一 :添加 模块 | 马 
Sine Wave Integrator Scope 


从 源 模块 库 〈Sources) 中 复制 正弦 波 模块 〈Sine Wave) 。 
连续 模块 库 (Continuous) 复制 积分 模块 〈JIntegrator ) 。 
输出 显示 模块 库 (Sinks) 复制 示波器 模块 〈Scope) 。 


步骤 二 ， 连 接 模块 总 四 加 


Sine Wave Integrator Scope 


步骤 三 ， 运 行 仿真 本 和 十 王 要 本 本 加 二 


-DLL 





四 


双击 示 波 带 模块 打开 Scope 窗 昕 双击 模 


型 窗口 菜单 中 的 【 Simulation>Start 】， 仿 和 











LT CT | 
让 TAO GO So 


如 果 将 以 上 算 例 的 初始 条 件 改 为 ，X(O) = 一 ] 


X(f) 一 SIn(L) 
x(O)= -1 入 组 初 如 条 作 C=0 


系统 的 解析 解 为 : “ X(1) 三 一 COS(T) 


Simulink 模 型 : 上 | :上 


Sine Wave Integrator Scope 
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单 目 由 度 系统 : 初始 条 件 : 


0 (OO 
要 求 ， 采用 Simulink 对 系统 进行 仿真 。 已 知 参数 :1 =]，c=1，A=1 
解析 解 为 : 


SIn(C 1) 


_ X 十 CCDX 
X(1)=e COS(C 四 十 as 


df 


人 
其 中 wo = 一 一 一 系统 固 有 频率 


mw, = mwML-5 ”一 一 阻尼 固有 频率 





en 一 一 一 一 相对 阻尼 系数 


11IX 十 CXY 十 KXY = EU 


Te/ 
Xe 从 已 知 参数 ;11 =] c=1 大 =1 
111 7171 


所 


c/m 


初 值 设置 为 0 初 值 设置 为 1 


于 二 站 寻 4 天 :下 


积分 得 到 速度 积分 得 到 位 移 
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如 果 系 统 中 没有 阻尼 ， 则 动力 方程 为 : 
11IX 十 KY = 已 知 参数 11 =1， 大 =1] 


初始 条 件 : 
人 


De 
解析 解 为 X(f) = X coS(Of 二 一 Sin(awo) 
WO 
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111 
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1.6 模型 基本 结构 
一 个 典型 的 Simulink 模型 包括 如 下 三 种 类 型 的 元 素 : 
图 信号 源 模块 
源 模块 系统 模块 显示 模块 
四 源 模块 广 -| 系统 模块 
国 输出 显示 模块 Simulink 模型 元 素 关 联 图 
> 信号 源 为 系统 的 输入 ， 它 包括 常数 信号 源 、 枯 数 信号 发 生 
器 〈 如 正弦 波 和 阶 跃 函数 波 等 ) 和 用 户 目 己 在 MATLAB 
中 创建 的 目 定 义 信号 。 
> “系统 模块 作为 中 心 模块 是 Simulink 仿真 建 模 所 要 解决 的 
主要 部 分 。 
> “系统 的 输出 由 显示 模块 接收 。 输 出 显示 的 形式 包括 图 形 显 
示 、 示 汲 器 显示 和 输出 到 文件 或 MATLAB 工作 空间 中 三 
种 。 输 出 模块 主要 在 Sinks 库 中 。 


SIMULINK (2) 


建 模 方 法 


本 章 内 容 和 学 习 目 的 


@ 本 章 介绍 一 些 模块 的 特殊 操作 ， 如 模块 操作 、 信 和 号 线 操 
作 、 模 型 注释 等 ， 对 模块 和 信和 号 线 有 一 个 整体 认识 ， 以 便 
为 后 续 的 动力 学 建 模 和 仿真 打下 基础 。 
@ 掌握 模块 操作 、 信 号 线 操 作 、 模 型 注释 的 基本 技巧 。 

2.1 打开 模型 

2.2 模块 操作 

2.3 信号 线 操作 

2.4 模型 注释 

2.5 创建 一 个 复杂 模型 

2.6 模型 和 图 形 捞 贝 到 Word 文 档 中 


-CH 下 
体 - 修 晤 


方法 1， 通过 【File>Open】 命 令 ; 
方法 2: 在 MATLAB 命令 窗 中 直接 键入 模型 名 。 


2.2 模块 操作 
2.2.1 调整 模块 大 小 

选中 模块 ， 模 块 四 角 出 现 小 方块 ， 然 后 按 住 鼠 标 拖 电 
2.2.2 旋转 模块 


选中 模块 ， 然 后 选择 药 单 命令 【Format>Rotate】， 模 
块 将 顺 时 针 方 癌 旋转 90 度 。 


2.2.3 模块 的 内 部 复制 
方法 1: 先 按 Ctrl 键 ， 再 单 击 模块 
方法 2: 碳 键 拖 上 折 模 块 
方法 3，【Edit>Copy】 和 【Edit>Paste 了】 
方法 4: 选中 模块 ，Ctrl+C 复 制 ，Ctrl 十 V 粘 贴 


2.2.4 删除 模块 
方法 1: 选中 模块 ， 按 Delete 键 
方法 2: 选中 模块 ， 然 后 从 模型 窗口 菜单 中 选择 【Edit>Clear 了】 
方法 3: 选中 模块 ， 然 后 选择 菜单 【Edit>Cut】， 删 全 筋 贴 板 
225 本 由 过 个 坑 专 
方法 1: 按 住 Shift 键 ， 同 时 用 雇 标 单 击 想 要 选中 的 模块 


方法 2: 使 用 “范围 框 ?。 在 使 用 “范围 框 ? 后 ， 者 像 继 续 选 择 
模块 ， 可 按 住 Shift 键 ， 再 单 击 要 选 的 模块 即 可 。 


各 多 个 模块 被 选中 ， 则 这 些 模 块 可 以 被 当 作 一 个 整体 进行 操 
全 人 这 全 向 


2.2.6 改变 模块 的 标签 

在 标签 的 位 置 上 双击 记 标 ， 则 模块 标签 进行 编辑 状态 。 
编辑 完 标签 后 ， 在 标签 外 的 任意 位 置 上 单 击 鼠标 ， 则 新 的 合 
2.2.7 改变 标签 位 置 

选择 菜单 命令 【Format>Flip Name】， 模 块 的 标签 将 发 
生 翻 转 。 者 原 标 签 位 置 在 模块 的 左边 ， 则 翻转 到 右边 。 寿 原 
标签 位 置 在 模块 的 下 方 ， 则 翻转 到 上 方 。 
2.2.8 隐藏 标签 

选择 菜单 命令 【Format>Hide Name】， 则 模块 的 标签 从 
屏幕 上 消失 。 此 时 【了 Hide Name】 变 成 为 【Show Name】， 选 
中 它 可 以 重新 显示 标签 。 
2.2.8 增加 阴影 

选择 菜单 命令 【Format>Show Drop Shadow】， 则 模 黄 
增加 阴影 。 此 时 【Show Drop Shadow】 变 成 为 【Hide Drop 
Shadow】， 选 中 它 可 以 隐藏 阴影 。 


2.3 信号 线 操 作 
2.3.1 绘制 信号 线 
2.3.2 移动 线段 
2.3.3 移动 节点 
2.3.4 删除 信号 线 
2.3.S 分 割 信号 线 置 
2.3.6 信号 线 标 签 


2.4 模型 注释 
2.4.1 增加 注释 的 方法 
2.4.2 改变 注释 字体 


注意 : 要 养 成 及 时 对 模块 和 信和 号 线 的 命名 以 及 及 时 对 模型 
加 标注 的 习惯 ， 这 会 提高 模型 的 可 读 性 。 


人 > 二 | | 《一 LA/ 亡 
光臣 和 的 咱 辣 下 全 


。 由 输出 端口 拖 电 鼠标 到 输入 端口 ， 或 拖 和 拐 鼠 标 由 输入 端 
口 到 输出 问 口 


Integrator Integrator 
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Gain Gain 


。 后 想 移动 信号 线 的 茶 段 ， 单 击 选 中 此 段 。 移 动 鼠 标 到 目标 
线段 上 ， 则 鼠标 的 形状 变 为 移动 图 标 。 按 住 鼠 标 ， 并 拖 鼻 
到 新 位 置 。 放 开 鼠 标 ， 则 信和 号 线 被 移动 到 新 的 位 置 。 


呈 可 
守 Integrator Integrator 
Et > CE 


Gain 






Gain 


。 知 想 移动 信号 线 的 茶 段 ， 单 击 选 中 此 段 。 移 动 鼠 标 到 目标 
线段 上 ， 则 鼠标 的 形状 变 为 移动 图 标 。 按 住 鼠 标 ， 并 拖 鼻 
到 新 位 置 。 放 开 鼠 标 ， 则 信和 号 线 被 移动 到 新 的 位 置 。 






Integrator 





Gain Gain 


移 J 册 | KLss 《- 
/ L 由 ?一 ] 可 
有 YY， 六 1 ee 2 g 


。 同 删除 模 坎 一样， 删除 信号 线 可 以 选中 信号 线 ， 然 后 按 


Delete 键 ， 或 者 利用 菜单 【EdiVyClear】 或 【EdiVyCut 】 
选项 进行 删除 。 


。 先 选 中 信号 线 。 按 住 Shift 键 ， 然 后 在 信和 号 线 上 需要 分 制 
的 点 上 单 击 也 标 ， 拖 动 节 所 到 新 的 位 置 ， 放 开 鼠 标 即 可 。 


加 > B 
Integrator Gain Integrator 


Gain 


@@) 


。 每 段 信 号 线 都 可 以 有 一 个 标签 。 双 击 要 标注 的 信号 线 ， 则 
言 号 线 的 附近 就 会 出 现 一 个 编辑 区 ， 在 编辑 区 内 输入 标签 
的 内 容 即 可 。 






Integrator 


Gain 


2.4 模型 注释 


使 用 模型 注释 可 以 使 模型 更 易 谍 懂 ， 其 效果 如 同 MATLAB 
程序 中 的 注释 行 一 样 。 对 于 经 常 使 用 Simulink 的 用 户 ， 养 
成 经 党 使 用 注 释 的 习惯 是 非常 重要 的 。 


2.4.1 增加 注释 的 方法 


在 模型 窗口 中 的 任何 想 要 注 酸 的 部 位 上 双击 鼠标 ， 将 会 出 现 
一 个 编辑 框 ， 在 该 框 内 输入 想 要 注释 的 内 容 即 可 。 
2.4.2 改变 注 酸 字体 


要 改变 注释 内 容 的 字体 ， 先 选中 注释 ， 选 择 模 型 窗口 沫 单 中 
的 【Format Fontj】 选 项 ， 就 会 出 现 一 个 字体 选择 的 对 话 框 
， 选 中 认为 合适 的 字体 ， 然 后 按 【OK]】。 


2.$ 创建 一 个 复杂 模型 
通过 两 个 算 例 来 说 明 建 模 中 的 其 他 一 些 技巧 


例子 1 


一 个 生长 在 缸 中 的 细菌 的 简单 模型 。 假 定 细菌 的 出 生 率 和 当前 细菌 的 总 数 成 
正比 ， 死 亡 率 和 当前 的 总 数 的 平方 成 正比 。 若 以 xz 代表 当前 细菌 的 总 数 ， 则 细 戎 


的 出 生 率 可 表示 为 : 
和 pit _ 7jale=DX 
细 兰 的 死亡 率 可 表示 为 : 
aedaln 7ale = DX 


细菌 总 数 的 总 变化 率 可 表示 为 出 生 率 与 死亡 率 之 差 。 因 此 系统 可 表示 为 如 下 
的 微分 方程 形式 ; 
X=DX 一 PX 


假定 六 =1/ 户 ,， 忆 =0.21/A ， 当 前 细菌 的 总 数 为 100， 计 算 一 个 小 时 后 镀 中 的 
细 苗 总 数 。 


元 = DX 一 DX- 


这 是 一 个 一 阶 系统 ， 因 此 一 个 用 来 解 微分 方程 的 积分 模块 是 必 
要 的 。 积 分 模块 的 输入 为 % 输出 为 x。 





和 X 
> = 一 > 
S 


Integrator 

















元 = DX 一 PX- 






























































增益 模块 ， 来 源 于 数学 乘法 模块 ， 来 源 于 数学 模块 
模块 库 〈Math ) 库 (Math) ， 用 于 实现 X 
与 
< 图 六 ] 
求 和 模块 ， 来 源 于 数 Gain Product 
学 模块 库 〈Math ) 
ww 
EN ] > 
Sum Integrator Scope 
《二 
示波器 模块 ， 来 源 于 输 
Wan 出 显示 模块 库 〈Sink) 


增 蔡 模块 ， 来 源 于 数学 积分 模块 ， 来 源 于 连续 
模块 库 〈Math ) 模块 库 〈Continuous ) 





































































于 1 
二 
Sum Integrator Scope 





Gain1 


。 久 

X=DX 一 PX 六 =1 PP=0.5 

当前 细菌 的 总 数 为 100， 计 算 一 个 消失 后 经 中 的 细菌 总 数 
步骤 三 : 设置 仿 其 参数 ， ii 


< 全 本 


双击 积分 模块 ， 
发 村 为 NE Product 将 初 信 设置 为 100 





















































Scope 





模型 的 起 始 时 间 默 认为 0， 终 止 时 间 默 认为 10.0。 需 要 改变 终止 
时 间 时 ， 旋 转 模 型 菜单 【Simulink>Parameters】， 打 开 模 型 参 
数 对 话 框 ， 设 置 【Stop time 】 为 10.0。 


现在 模型 就 全 部 完成 了 ， 选 择 【File>Save】 命 令 保存 模型 为 
bio_example，Simulink 将 以 bio_example.mdl 为 文件 名 保存 到 
指定 的 位 置 。 


| A LA 信人 一 个 一 | 
多 8 NT 一 = 


双击 示波器 模块 ， 并 选择 【Simulink>Start】 命 令 运 行 模 
汲 器 将 绘制 出 仿真 结果 。 





1 宣 


例子 2: 三 目 由 度 结构 的 强迫 振动 


世 计 请 果 “-: 





动力 学 方程 : 

1 0U 0012 3 一 上 0 1 书 SI o 〇 站) 

0 0 入 儿 交 0 -KK 3 |za 0 
局 SI O 〇 1) 

写成 矩阵 形式 : MI 不 +KXK = 五 () 环 (1) = 0 

0 

罗 3K -天 0 

村 才 本 人 作 下 作 全 全 全 人 

0 0 - 证 








MU 01 SEC OU 感 SIn(CC1) 
人 0 
0 0 7 玉 | 罗 0 一 3 2 0 


已 知 参 数 ， mm=1，/1, P=1,， 中 =0.5 


解析 解 : 倪 振 华 《 振 动力 学 》P189 P233 


xf (人 1 8 的 
OilenOUssli2 ssid jn2oa 0 全 sin(wD+021 = 全 sinw) 
ee 1 V3 国 





要 求 : 
采用 Simulink 对 系统 进行 动态 仿真 ， 并 与 解析 解 进行 对 
有 照 ， 计算 时 间 0 50 


局 SI o 〇 1) 
0 则 | 国人 0 
0 0 717|2X3 0 一 KK 38 | 5 0 
MIX +TKRX =P() 已 知 参数 : mm=1，/=1， 
X=-MTIRKRXE +HMTPO) P=1， 〇 = 0.5 






[0;0;0] 质量 阵 的 逆 M-1! 刚度 阵 K 


Scope 


质量 阵 的 逆 M- [1;0;0] Sine Wave ea 
频率 调整 为 0.5 
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2.6 模型 和 图 形 拷贝 到 Word 文 档 中 
2.6.1 拷贝 到 Word 文 档 


使 用 【Edit>Copy Model】 命 令 ， 然 后 在 Word 中 粘贴 


2.6.2 在 Word 文 档 中 编辑 模型 和 图 形 
使 用 取消 组 合 命令 


SIMULINK (3) 


运行 仿真 


介绍 两 种 Simulink 运行 仿真 的 方法 
3.1 使 用 窗口 运行 仿真 
3.2 使 用 MATLAB 命令 运行 仿真 


@ 和 营 握 以 上 两 种 运行 仿真 的 方法 


3.1 使 用 窗口 运行 仿真 


优 点 : 人 机 交互 性 强 ， 不 必 记 住 繁琐 的 命令 语句 即 可 进行 操 
作 。 使 用 窗口 运行 仿真 主要 可 以 完成 以 下 一 些 操作 。 


1. 设置 仿真 参数 ”仿真 参数 和 算法 选择 的 设置 


2. 应 用 仿真 参数 ”仿真 参数 和 算法 设置 后 ， 使 之 生效 


3. 局 动 仿真 选择 命令 运行 仿真 

4. 停止 仿真 选择 命令 停止 仿真 

S. 中 断 仿 真 可 以 在 中 断 点 继续 启动 仿真 ， 而 停止 仿真 则 不 能 
6. 仿真 诊断 在 仿真 中 若 出 现 错误 ，Simulink 将 会 终止 仿真 并 


在 仿真 诊断 对 话 框 中 显示 错误 信息 


选择 菜单 选项 【Simulation>Parameters】， 可 以 对 仿真 
参数 及 算法 进行 设置 ， 共 有 五 个 选项 卡 
。 解 法 设置 〈Solver ) 
。 工 作 间 IO (Workspace LO ) 
。 诊 断 页 (Diagnostics ) 
。 高 级 设置 (Advanced ) 
。 实 时 工具 对 话 框 〈( 了 Real-Time Workshop ) 


。 解 法 设置 (Solver) 〈 讲 ) 
设置 起 始 和 终止 时 间 ， 选 择 积分 分 解法 ， 指 定 求解 参数 和 
选择 输出 选项 
。 工 作 间 LIMO (WorkspaceLIJO) “〈 讲 ) 
管理 MATLAB 工 作 间 的 输入 输出 项 
。 诊 上 断 页 (Diagnostics) 〈 不 讲 ， 目 学 ) 
选择 在 仿真 中 警告 信息 的 等 级 
。 高 级 设置 (Advanced) (不 讲 ， 自 学 ) 
对 仿真 的 一 些 高 级 配置 进行 设置 
。 实 时 工具 对 话 框 〈( 了 Real-Time Workshop) (不 讲 ， 自 学 ) 
对 实时 工具 中 者 干 参数 进行 设置 。 奉 没有 安装 实时 工具 ， 不 
出 现 此 框 。 


当选 中 菜单 选项 【Simulation>Parameters】 后 ， 出 现 参 数 
及 算法 等 设置 页 。 再 点 击 【Solver】， 则 出 现 解法 设置 页 。 解 
法 设置 页 包括 三 项 内 容 : 


> 设置 仿真 的 月 动 时 间 和 终止 时 间 
> 选择 算法 并 指定 参数 
> 选择 输出 项 


仿真 时 间 


仿真 解法 〈 各 种 解法 说 明 见 下 页 ) 

默认 解法 〈ode45) 

变 步 长 解法 : ode4$S，ode23，odel13，odel1S，discrete 
定 步 长 解法 : ode$s，ode4，ode3，ode2，odel1，discrete 


最 大 步 长 
初始 步 长 


答 出 选项 : 用 户 用 来 控制 仿真 输出 个 数 的 对 话 框 ， 共 有 三 
个 菜单 选项 ， 定 义 输出 ， 产 生 附 加 输出 ， 产 生 
指定 输出 。 


各 种 ode 命 令 的 说 明 


型 
ode45 | 非 刚性 采用 4、5 阶 RungeKutta 法 | 大 多 数 场合 的 首选 算法 


非 刚 性 | 多 步 法 ， 采 用 Adams 算 法 ， 高 | ode45 计 算 时 间 太 长 时 
低 精 度 均 可 〈10-3 一 10-6) 取代 ode45 


[Ce 采用 梯形 法 则 算法 适度 刚性 


采用 2 阶 Rosenbrock | 当 ode45 失 败 时 使 用 ; 
度 中 等 或 存在 质量 矩阵 时 


适度 刚 和 
性 
算式 ， 精 
刚性 “| 一步 法 ， 采 用 2 阶 Rosenbrock | 低 精 度 时 ， 比 ode15s 
算式 ， 低 精度 人 或 存在 质量 和 矩阵 
两 阶段 算法 ， 低 精度 效 ; 或 存在 质量 矩阵 时 





(2) 工作 间 IO 


Simulink 作为 MATLAB 的 一 个 附件 ， 理 应 与 MATLAB 
很 好 地 结合 ， 它 的 输入 数据 可 以 从 MATLAB 的 工作 空间 中 获 
得 ， 其 仿真 结果 也 可 以 被 引入 到 MATLAB 的 工作 空间 。 实 现 
此 项 功能 需要 用 到 仿真 参数 对 话 框 中 的 工作 间 IO 页 。 


工作 间 IO 页 大 致 可 以 分 为 三 个 部 分 : 
分 从 MATLAB 工 作 空 间 获 得 系统 输入 〈Load from workspace) 
@@ 仿真 结果 输出 到 MATLAB 的 工作 空间 〈Save to workspace ) 
@ 输出 选项 〈Save option ) 


下 面 对 Workspace IO 选项 页 的 功能 与 使 用 分 别 予 以 简介 。 


争 从 MATLAB 工 作 空间 加 载 
(Load from workspace ) 
虽然 Simulink 提供 了 多 种 系统 输入 信和 号， 但 并 不 能 完全 满 
足 需 要 。Simulink 人 允许 使 用 用 户 自 定 义 的 信号 作为 系统 输入 信 
号 。 在 Load from workspace 框 中 ， 用 户 可 以 设置 MATLAB 
中 的 变量 作为 系统 输入 信和 号 或 系统 状态 初 什 ， 如 下 所 述 : 

(1) Input:， 用 来 设置 系统 输入 信和 号， 其 格式 为 [LuU， 其 中 
fu 均 为 列 回 量 ，ft 为 输入 信和 号 的 时 间 回 量 ，u 为 相应 时 刻 的 
言 号 取 值 。 可 以 使 用 多 个 信号 输入 ， 如 [6 ul，u2]。 输 入 信和 号 
与 Simulink 的 接口 由 Inport 模块 〈In1 模 块 ) 实现 。 


站， 
In1 Integrato Outi 口 


1 Scope 
[2 ER 上 2 
ms< Integrator1 Out 


扫 从 MATLAB 工 作 空 间 加 载 


(Load from workspace ) 


虽然 Simulink 提供 了 多 种 系统 输入 信号 ， 但 并 不 能 完全 满 
足 需 要 。Simulink 人 允许 使 用 用 户 自 定 义 的 信号 作为 系统 输入 信 
号 。 在 Load from workspace 框 中 ， 用 户 可 以 设置 MATLAB 
中 的 变量 作为 系统 输入 信和 号 或 系统 状态 初 值 ， 如 下 所 述 : 


(1) Input: 用 来 设置 系统 输入 信号 ， 其 格式 为 [5 ug， 其 中 
ft、u 均 为 列 回 量 ，t 为 输入 信号 的 时 间 回 量 ，u 为 相应 时 刻 的 
言 号 取 值 。 可 以 使 用 多 个 信号 输入 ， 如 [b ul，u2]。 输 入 信和 号 
与 Simulink 的 接口 由 Inport 模块 〈In1 模 块 ) 实现 。 

(2) JInitial state: 用 来 设置 系统 状态 变量 的 初始 值 。 初 始 值 
xnitial 可 为 列 问 量 。 
注意 : 使 用 xInitial state 所 设置 的 状态 变量 初始 值 会 目 动 履 
盖 系 统 模块 中 的 设置 。 另 外 ， 输 入 信号 与 状态 变量 需要 按照 系统 
移 理 中 Inport 模块 〈 即 Inl 模块 ) 的 顺序 进行 正确 设置 。 


急 仿真 结果 输出 到 MATLAB 的 工作 空间 


(Save to woOrkspace ) 

使 用 Workspace IO 选项 页 可 以 将 系统 的 仿真 结果 、 系 统 
仿真 时 刻 、 系 统 中 的 状态 或 指定 的 信号 输出 到 MATLAB 的 工 
作 空 间 中 ， 以 便 用 户 对 其 进行 定量 分 析 ， 如 下 所 述 : 

(1) Timne: 输出 系统 仿真 时 刻 。 

(2) States: 输出 系统 模型 中 的 所 有 状态 变量 。 

(3) Output: 输出 系统 模型 中 的 所 有 由 Output 模块 〈 即 
Outl 模块 ) 表示 的 信和 号 。 

(4) Final state:; 输出 系统 模型 中 的 最 终 状 态 变量 取 值 ， 即 最 
后 仿真 时 刻 处 的 状态 值 。 Ci | 他 


In1 Integrato Outi 


Scope 


In2 Integrator1 Out 


分 输出 选项 (Save option ) 


(1) Limit data points to last:， 表示 输出 数据 的 长 度 〈 从 信号 
的 最 后 数据 点 记 起 ) 。 


(2) Format:， 表 示 输 出 数据 类 型 。 共 有 三 种 形式 :Structure 
with Time 〈 带 有 仿真 时 间 变 量 的 结构 体 ) 、Structure 〈 不 带 
仿真 时 间 变 量 的 结构 体 ) 、Array (信号 数组 ) 。 


举例 说 明 : 


CT 疡 | 1 卜 C) 


In1 Integrato ee 口 
Scope 
到 
In2 = Out2 


Integrator1 


十 训 ， 
In1 Integrato outi 


2 1 | 


ms< Integrator1 


所 王 


Scope 


Out2 
假定 两 个 输入 信号 都 为 正弦 信号 sin(b， 并 且 假定 初始 值 为 [0，1] 
7 二 | sin( 人 = 一 cos(D+C 


yL = 一 cos( 雪 二 1] 
y2 = 三 一 COS( 世 十 2 


CD 一 站 卫 上 CT 


In1 Integrato 2 | 号 已 = 一 cos( 人 上 +1 
Scope y2 = 一 COS(f) 十 2 
2 | 
In2 Out2 


Integrator1 


将 Load from workspace 栏 的 两 项 全 部 选 上 ， 并 且 将 Input 栏 
改 与 为 [b u, U。 


将 Save to workspace 栏 的 四 项 全 部 选 上 。 
Save option 栏 的 三 项 分 别 为 : 1000，1，Array 


运行 仿真 前 ， 首先 需 要 生成 系统 短信 信 呈 与 状态 初始 值 ， 在 
MATLAB 命令 窗口 中 键入 如 下 命令 : 


>t=0: 0.1: 10; t= 旨 ; 


>U=Sn(t); 





> XInitial = [0， 了 如 ; 


GE 攻 和 一 


In1 Integrato E 加 y] 一 一 COS(T ) 平 | 
Scope y2 = 一 cosS(t) 十 2 
CD 一 | 革 全 一 >C2) 
In2 Out2 


Integrator1 





运行 whos 的 结果 : 


Name 91Ze Bytes Class 
f 101x1 808 double array 
tout 51X1l 408 double array 
U 101x1 808 double array 
XEFinal ID 10 double array 
XJnltlal < 10 double array 
XOUL 9 x2 810 double array 
yOUL | 816 doubjle array 


上 面 4 和 tout 的 维 数 不 相 同 ， 这 是 因为 在 Solver 中 采用 了 变 
步 长 解法 ， 大 采用 定 步 长 解法 〈 步 长 取 0.1) 则 维 数 相 同 。 


3.2 使 用 MATLAB 命 令 运 行 仿 芮 


MATLAB 提供 了 sim 命令 ， 用 户 可 以 在 MATLAB 的 环境 下 
以 命令 行 或 M 文件 的 形式 运行 Simulink 模型 。 


使 用 命令 行 方式 ， 用 户 可 以 在 脚本 文件 中 重复 地 对 同一 系统 在 
不 同 的 仿真 参数 或 不 同 的 系统 模块 参数 下 进行 仿真 ， 而 无 需 一 
次 又 一 次 启动 Simulink 图 形 窗 口中 的 Start Simulink。 


如 果 需 要 分 析 茶 一 参数 对 系统 仿真 结果 的 影响 ， 用 户 可 以 很 容 
易 地 通过 for 循 环 目 动 修改 任意 指定 的 参数 。 


下 面 对 进 行动 态 系统 仿真 的 命令 逐一 介绍 。 上 首先 看 一 个 例子 。 


建立 一 个 简单 的 动态 系统 ， 其 功能 如 下 : 
(1) 系统 的 输入 为 一 单位 幅 值 、 单 位 频率 的 正弦 信号; 
〈2) 系统 的 输出 信号 为 输入 信和 号 的 积分 。 


Sin(D) | sin(D)di 
要 求 如 下 : 
(1) 系统 的 输入 信号 由 MATLAB 工作 空间 中 的 变量 提供 ， 时 间 
0 一 10s; 


(2) 使 用 MATLAB 绘制 原始 输入 信和 号 与 系统 运算 绪 果 的 曲线 。 


问题 摘 述 : 


原 函 数 : Sin( 让 ) 时 间 区 域 : 1=0~ 10% 
原 函 数 的 积分 为 x(1) = | sin(Ddt = -cos(D+C 
若 取 零 初 始 条 件 ， 则 有 : ” C=1 


因此 ， 有 : xxX(1i) = 三 一 COS(I) 十 ] 


Xi) 三 一 CoOS(t 和 十 ] 


CL 全 
In1 Outi 
Integrator 
Simulink 模型 


(保存 文件 名 为 :command in _ out) 


对 Workspace IO 页 进行 如 下 设置 : 


(1) Load from workspace 栏 : Input 打 勾 ， 并 填 入 sim_input 
(由 MATLAB 工作 空间 输入 的 输入 名 ) 。Initial state 不 选 。 


(2) Save to workspace 栏 ， Time 和 Output 项 勾 上 ，State 和 
Final state 项 不 义 。 


(3) Save options 栏 依次 为 : 1000，1，Array 


X(1) 三 一 COS(1) 十 ] 


CT1) 人 CT) 
In1 Out1 
Integrator 


在 MATLAB 工作 空间 中 定义 输入 变量 sim_input 如 下 : 


>t= 0: 0.1: 10; t=tb; 匈 表示 输入 信号 的 时 间 范 围 
>U=Sin(t); % 产生 输入 正弦 信号 
>Sim_ input= [tu j]; % 传递 给 Simulink 系统 模型 的 变量 


然后 运行 Simulink 仿真 。 最 后 在 MATLAB 工作 空间 运行 如 下 


plot (ft u, tout, yout，-- ); 2rid 
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& 一 SIn(L) 


yOVULL = 三 一 COS(1oOV1) 十] 本 


人 











3.2.1 使 用 sim 命令 进行 动态 系统 仿真 
1 使 用 语法 
sim 命令 的 格式 为 : 


[tx y]=sim (model，timespan，options，nut ) 


[t，X， yl1，y2, ... ,yn ]= sim (model, timespan，options，nut ) 





以 上 是 完整 的 语法 格式 ， 实 际 使 用 时 可 以 省 略 其 中 的 某 些 参 
数 设 置 而 采用 默认 参数 。 除 了 参数 “mode， 外 ， 其 它 的 仿真 参数 
设置 均 可 以 取 值 为 空 抢 阵 ， 此 时 sim 命令 对 没有 设置 的 仿真 参数 
使 用 默认 的 参数 值 进行 仿真 ， 默 认 的 参数 值 由 系统 模型 框图 决定 
。 用 户 可 以 使 用 sim 命令 的 options 参数 对 可 选 参数 进行 设置 ， 这 
样 设置 的 仿真 参数 将 履 盖 模型 默认 的 参数 。 


2 参数 说 明 
sim 命令 的 格式 为 : 


[tx y]=simn (model，timespan，options，nutt ) 





[t，Xx，y1，y2, ..., yn ] = sim ( model, timespan，options，ut ) 


(1) model:， 需要 进行 仿真 的 系统 模型 框图 名 称 ; 


(2) timespan: 系统 仿真 的 时 间 范 围 《〈 起 始 至 终止 时 间 ) ， 可 有 
如 下 形式 : 


tFinal:， 设置 仿真 终止 时 间 。 仿 真 起 始 时 间 默 认为 0; 
[tStart tFinall]: 设置 起 始 时 间 (tsStarb 与 终止 时 间 〈(tFinal) ; 


[tStart OutputTimes tFinall]: 设置 起 始 时 间 (〈tStart) 与 终止 时 
间 〈tFinal) ， 并 且 设 置 仿真 返回 的 时 间 辐 量 [tStart OutputTimes 
tFinall]， 其 中 tStart、OutputTimes、tFinal 必须 按照 升序 排列 。 


sim 命令 的 格式 为 : 


[tx y]=simn (model，timespan，options，nutt ) 





[t，X， yl1，y2, ... ,yn ]= sim (model, timespan，options，nutt ) 


(3) options: 由 simset 命令 所 设置 的 除 仿 真 时 间 外 的 仿真 参数 ; 
(4) ut: 表示 系统 模型 项 层 的 外 部 可 选 输入 。ut 可 以 是 MATLAB 
函数 。 可 以 使 用 多 个 外 部 输入 utl、ut2、.…. 

(S$) ft: 返回 系统 仿真 的 时 间 癌 量 。 

(6) x: 返回 系统 仿真 的 状态 变量 矩阵 。 

(7) y: 返回 系统 仿真 的 输出 矩阵 。 按 照 顶层 输出 Outport 模块 的 
顺序 输出 。 如 果 输 出 信号 为 癌 量 输出 ， 则 输出 信号 具有 与 此 回 量 相 
同 的 维 数 。 


(8) y1，.…，yn: 返回 多 个 系统 仿真 的 输出 。 


> 乱 首 让 
3 举例 之 一 ， 人 简单 仿真 @ 一 四 一 全 
输入 : sin( 
X(1) 三 一 COS(1) 十 ] 文件 名 : _ command in out 


对 于 前 面 的 动态 系统 command_in_out， 在 前 面 进行 仿真 时 没 
有 使 用 命令 行 方式 ， 在 此 采用 命令 行 语 名 进行 仿真 。 在 仿真 之 前 上 
， 首 先 使 用 仿 0 话 框 设置 参数 ， 然 后 在 MATLAB 命 
令 窗 口中 键入 如 下 命令 : 


>t= 0: 0.1: 10; t=b; 4 表示 输入 信号 的 时 间 范 围 
>U=sin (ft); % 产生 输入 正弦 信号 
>sim_input=[bu];i 匈 传递 给 Simulink 系统 模型 的 变量 


>[ tout, x, yout ] = Sim (“command in _ out” ) 
% 使 用 sim 进行 系统 仿真 ， 仿 真 参数 取 与 前 面相 同 
>plot ( t u, tout, yout，“--” ); 2rid 





4 举例 之 二 : 仿真 时 间 设 置 


在 前 面 已 经 对 sim 命令 中 的 仿真 时 间 参 数 timespan 设置 
做 了 介绍 。 Timespan 具有 三 种 使 用 形式 ， 根 据 不 同 动态 系 统 
的 不 同 要 求 ， 用 户 可 以 选择 使 用 如 下 所 示 的 三 种 形式 进行 系 


统 仿真 : 
[txXy]|=sm (modeltkinal) 
[txy| =sm (model[ttart，tkinal] ) 





[t, xy =sm (model [tstart outputIimes tkinal ] ) 


[txXy]|=sm (modeltkinal) 
[txy]| =sm (model,[ttart，tkinal] ) 


[t xy =sm (model [tstart outputIimes tkinal ] ) 





此 外 在 默认 情况 下 ， 系 统 仿真 的 输出 结果 〈 输 出 时 间 、 状 态 和 
运算 结果 ) 受到 Simulink 求解 露 仿真 步 长 的 控制 ， 因 而 系统 仿 
真 输出 结果 也 受到 求解 露 步 长 的 控制 。 


如 果 需 要 在 指定 的 时 刻 输出 系统 仿真 结果 ， 则 需要 使 用 仿真 时 
间 设 置 的 第 三 种 方式 ， 其 中 [tsStart outputTimes tFinal] 表示 输 
出 时 间 同 量 ， 此 回 量 为 一 递增 行 问 量 。 在 sim 命令 中 使 用 
timespan 设置 系统 仿真 时 间 范 围 ， 会 覆盖 Simulink 原来 的 仿真 
时 间 设 置 ， 但 是 并 不 会 影响 到 系统 模型 。 


CT) 全 CCL 
In1 Out1 
Integrator 


输入 : sin(0) X(1) 三 一 COS(1) 十] 


仍 以 前 面 的 command_in_out 为 例 进 行 说 明 。 仿 真 参 数 设 置 对 
nn 使 用 四 组 不 同 的 仿真 时 间 对 此 系统 
进行 仿真 。 


为 了 与 系统 模型 文件 command_in_outmdl 区 别 ， 所 编制 的 M 
文件 名 取 为 ”command in out mm。 和 各 MX 文件 名 取 为 
command in_ outm， 则 在 MATLAB 工作 空间 运行 时 ， 将 返回 
系统 模型 文件 框图 。 


M 文 件 名 : command in_ out m 


>t=0:0.1:10; t=f;  u=sin(b); sim_ input=[tu]; 
>[toutlLxlyoutl]=sim(ccommand in_out?,S); % 仿真 时 间 范 围 0 一 Ss， 输 出 时 
间 回 量 toutl 由 Simulimk 的 求解 器 步 长 决定 
>[tout2,x2,yout2]=simn(ccommand in_ out',[1 8); % 仿真 时 间 范 围 1 一 8s， 输 
出 时 间 向 量 tout2 由 Simulimk 的 求解 器 步 长 决定 
>[tout3,x3,yout3]=sim(ccommand in_out',[1:8); % 仿真 时 间 范 围 1 一 8s， 并 
且 每 隔 ls 输出 一 次 ， 即 输出 时 间 变 量 为 [L234S678] 


>[tout4,x4,yout4]=simn(ccommand in_out?,1:0.2:8); % 仿真 时 间 范 围 1 一 8s， 
并 且 每 隔 0.2s 输出 一 次 ， 即 输出 时 间 变 量 为 [L1.2 1.4 .… 7.6 7.8 8] 


>Subplot(2,2,H); plot(tu,toutl,yout1, 关 7); 


>Subplot(2,2,2); plot(tu,tout2,yout2,”*#7); 
>Subplot(2,2,3); plot(tu,tout3,yout3, 衬 7); 





>Subplot(2,2,4); plot(tu,tout4,yout4,?*7); 


>t=0:0.1:10; t=tf; u=sin(D; sim_input=[bu]; 


>[toutl,xl,youtl=sim('command in _ out,S); 
>[tout2,x2,yout2]=sim( command in _ out,[1 8]); 


>[tout3,x3,yout3]=sim( command in _ out,[1:8]); 
>[tout4,x4,yout4]=sim(command_in_ ouft” ;1:0.2:8); 


>Subplot(2,2,1); plot(tu,toutlLyout1 ,> 关 7); 
>Subplot(2,2,2); plot(tu,tout2,yout2," 关 7); 





























































































































>Subplot(2,2,3); plot(tu,tout3,yout3,”*#7); 




























































































>Subplot(2,2,4); plot(tu,tout4,yout4,7); 
输入 : sin() 1 




















X(1) 三 一 COS(t) 十] 
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S$ 举例 之 三 : 外 部 输入 变量 设置 

前 面 对 动 态 系 统 command in_ out 进行 仿真 时 ， 通 过 设置 
Simulink 仿真 参数 设置 对 话 框 中 Workspace IO 中 的 外 部 变量 
和 输入， 以 使 系统 在 仿真 过 程 中 获取 输入 信号 sim_input。 除 了 使 
用 这 种 方法 从 MATLAB 工作 空间 中 获得 系统 输入 信号 之 外 ， 
用 户 还 可 以 通过 使 用 sim 命令 中 的 ut 参数 来 设置 系统 的 外 部 和 输 
入 信号 。 下 面 介 绍 如何 使 用 ut 参数 设置 外 部 输入 信和 号。 


仍 以 command in_out 为 例 进 行 说 明 。 使 用 如 下 命令 : 
[t x, y]=sim (model timespan, options, ut) 


ut 为 一 个 具有 两 列 的 矩阵 ， 第 一 列表 示 外 部 输入 信号 的 时 间 ， 
第 二 列 代表 与 时 间 列 相对 应 的 外 部 输入 信和 号 。 


输入 : sin( ES 


X(t) 三 一 COS(D) 十 ] | Integrator 2 


MATLAB 程序 〈 程 序 名 : command in_out_m2.m) 


>t=0:0.1:10;  t=tf; au=sin(D0; sim_input=[t，U]; 


>[tout1，x1，youtl]=sim(ccommand in_ out?，10); % 使 用 Simulink 仿真 参 
数 对 话 框 中 的 workspace IO 从 MATLAB 工作 空间 中 获得 输入 信和 号 


>U=coSs(t); 


>ut=[t U]; % 改变 系统 输入 信和 号 


>>[tout2，x2，yout2]=sim(ccommand_in_out?，10,，[]，ub; % 使 用 Sim 中 的 
ut 参数 获得 系统 输入 信号， ut 的 使 用 会 履 羡 由 Workspace LO 的 系统 输入 设 
置 ， 这 一 点 可 以 在 下 面 的 系统 仿真 结果 图 形 中 反映 出 来 


>Subplot(I，2， 了 也 ; “plot(tout1，youtl); ”crid 


>Subplot(I，2，2); “plot(tout2，yout2); “grid 





yoxll = | sin(toxtDd(oxtD = -cos(ioxrD+1 


yoxt2 = | cos(ioxt2)Jd(iox1i2) = sin(iox12]) 














人 
1.8 In1 Outi 
Integrator 
0.5 
yo 1 = | sin(ioxu DJd(iowdD = -cos(1oxTD+1 















































0 yo 2 = | cos(1ol 2)d(io 2) = sin( io 2) 









































图 中 上 图 表示 系统 输入 为 sin(O 时 的 相应 曲线 ， 下 图 表示 cos(b 时 的 相 
应 曲线 。 从 图 中 可 明显 看 出 ， 当 使 用 sim 命令 的 ut 参数 时 ，Simulink 仿真 
参数 设置 对 话 框 中 的 设置 被 覆盖 。 以 前 对 话 框 中 的 外 部 输入 是 名 为 
sim_input 的 正弦 信号 ， 而 采用 ut 参数 后 执行 的 余弦 输入 信号 。 注 意 : 这 里 
指 的 “覆盖 ”， 并 不 是 在 Workspace IO 对 话 框 的 Input 中 ， 将 sim_input 改 
0 ut， 事 实 上 并 没有 改变 ， 只 是 不 执行 sim_input， 而 执行 了 命令 行 中 
J ut。 


3.2.2 simset 命令 


simset 命令 是 用 来 创建 和 编辑 options 结构 的 。 


3.2.3 simget 命令 
simget 命令 用 来 获得 系统 模型 的 仿真 参数 设 


3.2.4 simplot 命令 
使 用 语法 : 
simplot ( data ); 
simplot ( time，data ); 


说 明 : data 是 动态 系统 方针 结果 的 输出 数据 ， 一 般 由 outport 模 
块 、To Workspace 模块 等 产生 的 输出 ， 其 数据 类 型 可 以 为 矩阵 
、 问 量 或 是 结构 体 等 。 


SIMULINK (4) 


基本 模块 介绍 


@ Simulink 模块 的 基本 知识 
@ 基本 模块 的 性 质 分 关 和 基本 操作 


4.1 连续 系统 模块 


大 多 数 物理 系统 可 以 用 微分 方程 进行 撒 述 ， 因 此 可 以 用 连 
续 系 统 模拟 。 最 简单 的 模型 是 线性 模型 和 定 第 模型 。 


例如 ， 振 动 理论 中 的 动力 学 方程 : 
hCGI)+TCX(I) 二 KZ = 书 (1 


其 中 ，x 为 系统 的 广义 坐标 列 巾 量 ，M 为 质量 窍 阵 ，C 为 阻 
尼 和 矩阵 ,天 为 刚度 年 阵 ，P(O 为 外 部 激励 列 向 量 。 


在 Simulink 中 ， 用 来 模拟 连续 系统 的 基本 模块 有 四 个 : 增 义 模 
块 ， 求 和 模块 ， 微 分 模块 ， 积 分 模块 。 除 了 这 四 个 基本 模块， 
传递 函数 模块 也 经 向 用 来 模拟 物理 系统 和 控制 需 。 


4 二、 不 于 二 岂 


作用 : 使 增益 模块 的 输入 信和 号 乘 以 一 个 常数 ， 并 输出 。 
可 用 代数 表达 式 表示 为 : 


y(D = 大 xD 
简 图 如 下 : 
1 yG 一 > 一 
Gain 
Simulink 增 义 模块 图 


增益 模块 
注意 : yD、x(O、 天 可 以 为 标量 、 辐 量 或 惩 阵 。 
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加 量 和 标量 的 乘积 ED 
[2| 运 2 瑟 [24| 安 EU 到 


Constant 二 Display 
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[ 2 专 [7 10| Constant Display 
Gain 
矩阵 和 回 量 的 乘积 1 
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Display 


Gain 


2. 求 和 模块 
作用 : 对 两 个 或 多 个 信号 进行 求 和 运算 。 
可 用 代数 表达 式 表示 为 : 


C=w+D 


两 种 形状 : 圆 形 和 方形 。 


0 C 于 
D 
求 和 模块 必须 至 少 有 一 个 输入 而 仅 有 一 个 输出 。 输 入 的 正 负 
号 的 数目 由 双击 模块 进入 编辑 栏 进行 设 定 。 


求 和 模块 个 但 可 以 进行 标量 求 和 和 运算， 也 可 以 进行 四 量 或 矩 
阵 求 和 运算 ， 但 是 标量 或 矩阵 的 维 数 必 须 相 等 。 


标量 求 和 Eu 
Constant D 
isplay 


易 人 3 一 | | s。 | constanti 


四 加 日 Constant 
下 Constanti 


窍 阵 求 和 1 3 


is | 十 
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Ai AN 上 此 了 二 


作用 : 计算 输入 对 时 间 的 变化 率 。 
axX 


代表 如 下 微分 方程 : 本 
微分 模块 如 图 所 示 : tn 


Derivative 


例如 ， 对 于 动力 学 方程 : 
LIXC1)+TCXGT+TKXCI)= 忆 (1) 


考虑 对 正 驴 信号 si (0 的 微分 ， 人 Sm = cos( 


Simulink 模型 框图 和 仿真 结果 如 下 页 图 形 所 示 。 


df[simn(GD)] 
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Sine Wave Derivative 














os0) | 引 


Scope 


从 _ -HH | 汪汪 
人 全 基 
刀 4 天 坟 


作用 : 计算 输入 信号 从 起 始 时 间 到 当前 时 刻 对 时 间 的 积分 。 
代表 如 下 微分 方程 y(D)= ydo)+| x(z)dz 








积分 模块 如 图 所 示 : 

例如 : Integrator 

y(D= | sin(Ddi 村 交 更 站 要 本 要 哆 罗 硕 

0 人 

nl 和 
四 国 了 硬 泛泛 二 本 和 省 荐 因 
二 Scope 0 








Sine Wave “Integrator 


初 值 设 置 为 -1 


S$. 简单 物理 模型 


利用 前 面 所 介绍 的 这 些 模块 可 以 模拟 由 线性 微分 工程 描述 的 任何 物 
理 模型 。 例 如 ， 考 虑 如 下 所 示 的 简单 的 小 车 系统 运动 。 
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有 忽略 摩 控 力 ， 过 Y 一 一 一 
0 





雯 





ER 加 入 两 个 积分 模块 ， 第 一 个 模块 用 来 计 
可 用 模块 图 表示 为 : 算 速度 ， 第 二 个 模块 用 来 计算 位 移 : 


人 区 丁目 
假定 F=sin() 为 正弦 激励 


，m=0.5。 求 0 一 20s 区 间 忆 关 四 
内 的 系统 位 移 响应 曲线 。 : 于: 册 








Sine Wave 











40 


20 


10 


若 要 求 同 时 输出 位 移 和 速度 ， 则 模型 框图 为 : 


Sine Wave 
Scope 

















吉 要 求 同 时 输出 位 移 、 速 度 和 加 速度 ， 则 模型 框图 为 : 





Sine Wave 

















11IX 十 CY 十 AX=U0，xXU)=0.2， XU) = 
1 一 0.9，K=1 ec=0.04 


CS 


OCRET 4 U.0 . 1 
X YX= 一 一 -一 X 一 
111 111 0.9 0.9 









Integrator Integrator1 ” Scope 人 | 


初 值 设置 为 0.2 


Gainl 人 











6. 传递 函数 模块 

传递 函数 表示 法 频 索 地 应 用 于 控制 系统 设计 和 系统 的 动态 模拟 。 传 
递 函 数 定义 为 系统 在 零 初 始 条 件 下 输出 的 Laplace 变换 与 输入 的 Laplace 
变换 之 比 。 传 递 函 数 是 一 种 描述 系统 动力 学 输入 输出 关系 的 简便 方法 。 








例如 : 
考虑 弹 筑 质 量 阻尼 系统 动力 学 方程 可 写 为 : 
1 人 学 十 CX 十 = 下 
“9 作 Iaplace 变换 ， 并 忽略 初始 条 件 ， 有 : 
nn 113“X(S)+CcSX(CJ+KXC) = 开 (y) 
则 系统 的 传递 函数 可 表示 为 : 
从 (3) ] 1/7 
wa 
(8) 118 十 CS 十 ai 


1717 1717 


111 十 CXY 十 Kx 三 广 


aa 
忆 (8S) 118 二 CS 二 + asia 
1 711 
使 用 基本 模块 元 素 所 建立 起 的 系统 模型 
相对 应 的 传递 函数 模型 


1/m 
s 4(c/m).Ss+(k/m) 


Transfer Fcn 





以 前 面 受 正弦 外 力作 
用 的 质量 块 为 例 :  .， 下 





111 
人 F=sin(b) 
Im=(0.S 
传递 函数 仿真 框图 : 
P > 3 
有 Scope 


Transfer Fcn 


传递 函数 分 子 输 入 : 10.5 
分 母 输入 : [1 0, 0] 


必 必 


Sine Wave 


| | | | 
0 2 4 6 8 


Simulink 模型 仿真 框图 为 : 





7 


人 


0 














对 传递 函数 模块 的 说 明 : 

“ 用 传递 函数 模块 对 线性 定 利 系统 进行 仿真 ， 能 够 使 仿真 模型 
简单 和 紧凑 ， 但 是 无 法 和 输出 内 部 变量 ， 如 x 的 导数 。 

“无 法 适用 具有 轧 始 条 件 的 情况 。 


“传递 郴 数 模 块 只 适用 于 单 输入 单 输出 系统 ， 即 单 自 由 度 系 统 
， 但 无 法 应 用 于 多 输入 多 输出 系统 ， 即 多 自由 度 系统 。 


4.2 状态 空间 模块 


状态 空间 模块 可 以 起 到 与 传递 函数 模块 相同 的 作用 。 上 所 不 同 
的 是 ， 状 态 守 间 模块 允许 用 户 指 定 初 始 条 件 ， 并 且 可 以 共 这 
着 


1. 状态 空间 的 概念 


考虑 如 下 单 自由 度 系统 :11IX(C1) 十 CX() 十 XI) = 尼 () 


xD 为 系统 的 三 义 坐 标 ; 大、 及 、c 分 别 为 单 自由 度 系统 的 质量 、 刚 度 、 阻 
尼 系 数 ; 下 为 系统 的 外 部 输入 。 这 些 变量 都 是 标量 。 


X(L) 
令 : 了 (站 = 0 以 X(t) 和 X(1) 构 成 的 空间 即 为 状态 空间 
7(D 为 一 个 2 维 列 同 量 ， 称 为 系统 的 状态 辐 量 。 利 用 状态 同 量 ZYOD， 则 上 述 
动力 方程 可 转化 为 : 


国 0 1 | 的 0 
Ca 有 ai 
X(T) --111 KE 一 1 cx 17 


(1) -4 (1) “和 痊 入 惩 阵 


系统 官 阵 


14X(D+CX (CD+ 开 和 (人 = 互 (人 
素 ( 人 =[X (恩人 (人 (ER 尺 “ 寺 维 广义 坐标 列 向 量 


h e 有 ”质量 矩阵 C < 尺 "入 阻尼 矩阵 
及 <R 民 5 刚 度 和 矩阵 不 (1 站)E 民 外 部 输入 列 矩 阵 


厌 ( 
令 : 了 (四 三 | - 岂 as 27 维系 统 状 态 列 回 量 
不 (人 
系统 的 状态 方程 : 7( 站 = 4 了 (站 上 + 妇 已 (CD) 
L/ 7 L/ 
4 芝 二 E Re 万 到 E se 
-MI - 开 -MC 1 


系统 官 阵 输入 矩阵 


2. 线性 单 输入 输出 系统 


炒 





1X 上 cz+K=0，x(0)=0.2，XFO)=0 


1 全 5 有 起 05 





C 





SN 





写成 状态 方程 形式 ， 有 :了 (1) = 4 了 (人 人) 


和 二 了 


在 状态 空间 模块 中 输入 : 

A=[0，1; -2，-0.1] 本 
B=[0; 2] %% 因 为 没有 外 部 输入 ， 也 可 B=[0; 0] 中 y=cem 
C=[L 1 % 不 能 写成 C=[0，0]， 否 则 输出 为 0 一 SEE 
D=0 


Scope 








上 二 一 一 一 一 一 一 一 咎 涪 上 一 一 一 一 一 一 一 一 一 一 圭一 一 一 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 一 了 一 一 一 一 一 一 一 一 一 一 


人 


和 


本 一 一 可 二 
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说 明 : 对 于 本 算 例 ， 因 为 考虑 

的 是 没有 外 部 输入 的 目 由 振动 

， 因 此 状态 空间 模块 可 以 没有 
Scope 左 端 输 入 。 


的 结果 。 各 取 C=[1,1]， 则 
在 Scope 的 图 形 中 ， 将 显 
示 出 位 移 和 速度 数据 全 加 
后 的 数据 的 图 形 。 


示波器 输出 为 y 的 图 形 


入 
和 


香 取 C=[1，0]， 则 示波器 显示 


呈 出 位 移 时 程 


Scope 0.15 








上 十 二 二 周三 二 二 二 站 区 呈 二 臣 二 区 居 二 匡 避 和 二 二 辣 熙 司 攻 训 二 二 二 品 


State-Space 











也 可 : C=[1.0;0,1H] ， D=zeros(2,T 
直接 得 出 以 下 图 形 奋 取 C=[0，IH， 则 示波器 显示 








Y 程 
YY 
0.3 速度 日 个 甩 
0.3 
本 | 
本 一 
四 
本 避 | 
人 
人 到 
四 
-AL----L---- 一 
四 
媚 守 1 工 _ | 
Ed | 
0.4 0.4 




















3. 多 输入 输出 系统 

11 三 U.3， 1 =] 
[| 
Ce 








SSE SN 


2 (1) 
MG 17 0 扩 十 人 CE C = Cl 十 C>? 一 (2 
0 177> 人 人 一 (2 人 


RN 四 本 | 日 
X(0) 0 X(0)| |0 


14XZCD+CXECD+RKRKE(CD=O 从 (人 ) nn 


11 证 王八 777 汪 二 中 
和 
可 二 0 由 


人 
1 语 (D)+C 写 (DTKKC)=O 三 na | 由 | 





JR 大和 双 一 4Y 1 了 OO 十 从 (1) 沽 这 / 72 
写成 状态 方程 : (1) 人 ( ) 书 ( 四 _ 1 天 _11RmC 


初始 条 件 : 系统 矩阵: 
0 0 0 1 0 
人 (0O) | | 0 0 0 0 1 
| 
RS 0 全 -6 4 -02 0 


0 ah 


0 0 1 0 0.2 


1IX = 
(+TCX(CDT+TKECD=O | | 0 
汪 7(0) = 
2 -2 0.05 -0.05 0 


在 状态 空间 模块 中 输入 : A=[0,0,1.0;0,0,0,1;-6,4,-0.2,0.1;2,-2,0.0S,-0.0S] 
B=[0;0;0;0] % 因 为 没有 外 部 输入 
C=[1.11.1] % 不 能 写成 C=[0,0,0, 0], 否 则 输出 为 0 
D=0 





Scope 


State-Space 





Scope 


State-Space 


在 MATLAB 工 作 宰 间 画 出 


示波器 显示 的 第 一 的 两 个 质量 的 位 移 时 程 
个 质量 的 位 移 时 程 也 可 : C=[L0,0,0;0,1.0,0] 








C=[1.0,0.0] D=zeros(4,1) 直 接 得 出 





到 


人 


二 二 河和 二 全 二 半生 二 二 扼 相生 敌 和 二 下 光线 条 是 二 二 二 折 冯 渗 
| | 


至 刺 苛 训 抽 知礼 党 有 


一 一 一 二 一 上 村 一 二 二 二 和 学 二 二 本地 二 本 二 二 二 二 王 富 近 三 一 全 二 二 记 - 林 一 呈 ' 安 二 二 ,二 ,一 


党 全 忆 后 三 ( 二 司 吉 计 王 呈 司 呈 和 中 已 


-rr----=---------+-------r--------------- -一 -一 一 一 一 


Scope 





State-Space 


在 MATLAB 工 作 空间 画 出 的 两 个 质量 的 速度 时 程 
也 可 : C=[0.0,1.0;0,;0.0,1]， D=zeros(4, 了 T) 直 接 得 出 


0.4 





和 











4.3 离散 系统 模块 


离散 系统 和 采用 差分 方程 表示 。 离 艇 信和 号 是 一 组 以 来 样 周期 为 
间隔 的 离散 时 间 序 列 。 对 于 大 多 数 的 物理 系统 ， 信 和 号 原本 都 是 在 
时 间 上 连续 的 。 但 在 对 信号 的 采集 过 程 中 ， 需 要 通过 传感器 等 采 
集 工 具 对 信和 号 进行 采集 ， 因 此 最 终 得 到 的 信号 都 是 以 采样 周期 为 
时 间 间 隅 的 离散 数据 。 寿 原 物 理 系 统 也 用 离散 数据 形式 进行 描述 
， 则 构成 了 离 若 时 间 系 统 。 


. 离散 增益 模块 
上 。 与 连续 系统 采用 的 模块 相同 


. 离散 求 和 模块 
. 离散 延迟 模块 
. 时 间 离 敢 积 分 模 其 
简单 离 敢 系统 模型 
. 离散 传递 函数 模块 


人 上 ww 王 


4.4 离散 状态 空间 模块 
1. 连续 时 间 系 统 的 离散 化 
考虑 慰 自由 度 线性 定常 系统 :了 (CD) = 47(D+Bu(D) 
革 示 ER Recee ie Re 
式 〈1) 的 解析 解 了 (CD =e 7do)+| e40-DBu(r)dz 
is al 和 其 中 : 7 为 数据 采样 周期 
并 采用 如 下 零 阶 保持 器 :CDD =xK)，H<r<(E+DT7 


(KE+JT) 


则 式 (3) 转化 为 ， (KE+D = e47 了 ( 提 + | e4IeDT7-ndtx( 有 


KT 


令 ， 1 三 (KE 二 JU7Z 一 7 
则 上 式 变 为 了 (KE+D=ee7(o+| ce4rd1Bu( 有 ) 


(1 ) 


〈2 1) 


《3 ) 


〈4) 


(CD) = 47(D)+ 万 xD) (1) 


7G+D=es(GD+| ce4d1Bue(A (4) 
令 ， 下 =e47 < Room， G=| edBesR” 
则 式 〈4) 转化 为 :了 (EC+D = FF(O+Gu(1) (6) 


式 〈$) 中 的 两 个 变量 将 于 有 限时 间 内 趋 于 各 值 矩 阵 。 


这 样 ， 式 〈1) 所 示 的 连续 时 间 系 统 就 转换 成 了 式 〈6) 所 示 的 离散 时 间 系 统 
。 怎 阵 正 称 为 离散 时 间 系 统 的 系统 算 阵 ，G 称 为 离散 时 间 系 统 的 输入 和 矩 阵 。 


2. 忆 和 6G 的 计算 
根据 矩阵 指数 的 定义 ， 有 ， 下 =e7 =14+4T+ 人 = 袜 人 


=0 











下 o0 0 oo IT 
同时 不 难 求 得 ，GCG= | edB = 已 = 万 
上 之 (上 十 ])! 上 人 


Xx(n+1)=Ax(n)+Bu(n) 
y(n)=CXx(n)+Du(n) 





Discrete State-Space 
4. 采用 离散 状态 空间 模块 进行 仿真 


对 于 两 目 由 度 系统 : 
MX(CD)+CXK(CD+KXCGD)=O 
和 (0)=[0.2 0]” X(OO)=[0 0 


2 卫 风 7 _[mm 0 En 
(人 0 1772 







14XCODHTCKCD+KE(CD=O X(O=[02 0 过 (0D)=[0 0 


> 站 太子 ; 入 ( 0 7 
居 候 下 到 (四 = . 和 2x2 2x2 
写成 状态 方程 ， 了 (0D) = 47(D) 了 CO) 并 4 _ Eee 国 


大 (0) 


初始 条 件 ，YO-| xzo| -2 oo 可 
离散 时 间 系统 的 系统 矩阵 〈 候 


连续 时 间 系 统 的 定 采 样 周期 为 了 7 王 0.02) 
系统 窍 阵 : (程序 名 : chapS_sectionS3.m) 
0 0 ] 0 0.9988 ”0.0008 “0.02 0 
is 0 0 0 ] 到 0.0004 “0.9990 0 0.02 
stat 人 RD 0.1 一 0.1197 ”0.0798 ”0.9948 0.0028 
) _》 005 005 0.0399 -0.0399 0.0014 0.9986 


(KK 二 TD = 正字 (A) 7(0)=[0.2 0 0 (0 


0.9988 ”0.0008 ”0.02 0 


yY(K+H) = 天 7Y(C) 0.0004 ”0.9990 0 0.02 
| -0.1197 0.0798 ”0.9948 0.0028 
7(0)=[0.2 0 0 0 


0.0399 -0.0399 0.0014 0.9986 
在 状态 空间 模块 中 输入 : A=[0.9988,0.0008,0.02,0;0.0004,0.9996,0,0.02;- 
0.1197,0.0798,0.9948,0.0028;0.0399,-0.0399,0.0014,0.99865] 

B=[0;0;0;0] % 因 为 没有 外 部 输入 

C=[111.1 % 不 能 写成 C=[0,0,0,0], 和 否则 输出 为 0 

D=0 
初 值 输入 为 : [0.2; 0; 0; 0] 
步 长 输入 为 002 Xx(n+1)=Ax(n)+Bu(n) 


y(n)=CXx(n)+Du(n) 


Scope 





Discrete State-Space 





示波器 显示 的 第 一 
个 质量 的 位 移 时 程 





| 


和 本 二 全 河 导 要 二 瑟 呈 二 二 机 号 司 有 二 本 二 二 尼 二 瑟 二 全 二 后 汪汪 二 后 辣 汪 二 二 


局 己 后 忆 所 局 了 局 后 局 后 二 局 剖 司 局 后 汪 后汉 志 二 二 二 语 











X(n+1)=Ax(n)+Bu(n) 


y(n)=CXx(n)+Du(n) 


Scope 


Discrete State-Space 


在 MATLAB 工 作 空 间 画 出 
的 两 个 质量 的 位 移 时 程 





上 = 


加 世 芝 本 必 必 人 


下 


HH--------------+-------r---------------+--- 一 一 -一 











与 前 面 连续 时 间 系 统 的 仿真 结果 一 致 。 


本 章 学 习 内 容 和 目的 


@ 进一步 掌握 连续 系统 的 建 模 和 仿真 技巧 
@ 竺 握 癌 星 线性 系统 的 描述 方法 
e 和 


于 


积分 模块 
高 级 积分 器 


在 使 用 Simulink 对 实际 的 动态 系统 进行 仿真 时 ， 积 分 
运算 可 以 说 是 Simulink 求解 器 的 核心 技术 之 一 。 在 前 面 章 
节 中 使 用 过 积分 模块 对 数据 进行 积分 处 理 ， 事 实 上 我 们 所 使 
用 的 是 积分 模块 的 一 种 简单 积分 方式 。 积 分 模块 还 有 一 种 积 
分 方式 一 一 重 置 积 分 方式 。 所 谓 的 重 置 积 分 方式 是 指 当 重 置 
言 号 触发 时 将 模块 的 输出 重 置 为 初始 条 件 。 本 节 积 分 模块 的 
重 置 积分 方式 一 一 高 级 积分 器 进行 介绍 。 








首先 对 积分 器 的 各 个 跨 口 进行 简单 的 介绍 。 


下 图 分 别 为 使 用 缺 省 参数 设置 下 的 积分 磊 外 观 与 选择 所 用 参 
数 设 置 下 的 积分 器 的 外 观 比较 。 





Integrator 


Integrator 





【External reset】 为 外 部 重 置 设 置 。 它 用 在 当 重 置信 号 发 生 触发 
事件 时 ， 模 块 将 按照 初始 条 件 重 置 状 态 量 。 可 以 采用 不 同 的 触发 方 
式 对 积分 器 状态 进行 重 置 : 


> none: 关闭 外 部 重 置 ; 
> rising: 当 模 块 接收 到 的 触发 信号 上 升 通过 零点 时 ， 重 置 过 程 开 


始 ; 
> falling: 当 模 块 接收 到 的 触发 信 呈 下降 通过 零点 时 ， 重 置 过 程 开 
始 ; 


> either: 无 论 触 发 信号 上 升 或 下 降 通 过 零点 ， 重 置 过 程 都 开始 ; 
> level: 当 触 发 信号 非 零 时 ， 使 得 积分 器 输出 保持 在 初始 状态 。 





【initial condition source】 为 初始 条 件 设 置 。 设 置 积分 
器 初始 条 件 的 方法 有 两 种 : 

> external: 从 外 部 输入 源 设 置 初 始 条 件 。 初 始 条 件 设置 

端口 以 x0 作为 标志 ; 

> internal: 在 积分 器 模块 参数 对 话 框 中 设置 初始 条 件 ， 说 


明 模 块 的 初始 值 是 从 内 部 获得 的 。 选 上 后 ， 下 面 将 出 现 要 
求 输入 初始 值 的 输入 栏 。Internal 为 默认 设置 。 





【Limit output】 积 分 器 输出 范围 限制 。 在 某 些 情 况 下 ， 积 分 器 的 
输出 可 能 会 超出 系统 本 身 所 允许 的 上 限 或 下 限 值 ， 选 择 积分 器 
输出 范围 限制 框 〈Limit output) 并 设置 上 限 值 CUpper 
saturation limit) 与 下 限 值 (Lower saturation limit) ， 可 
以 将 积分 器 的 输出 限制 在 一 个 给 定 的 范围 之 内 。 此 时 积分 器 的 
输出 服从 下 面 的 规则 ; 


> 当 积 分 结果 小 于 或 等 于 下 限 值 并 且 输 入 信和 号 为 负 ， 积 分 器 的 输 
出 保持 在 下 限 值 (下 饱和 区 ) : 

> 当 积 分 结果 在 上 限 值 与 下 限 值 之 间 时 ， 积 分 器 得 出 为 实际 积分 
结 傈 ; 


> 当 积 分 结果 大 于 或 等 于 上 限 值 并 且 输 入 信号 为 正 ， 积 分 器 的 输 
出 保持 在 上 限 值 (上 饱和 区 ) 。 





【Show saturation port】 在 积分 器 中 显示 饱和 端口 。 此 
端口 位 于 输出 端口 的 下 方 。 饱 和 端口 的 输出 有 三 种 情况 ， 用 
来 表示 积分 器 的 饱和 状态 : 


> 输出 为 1， 表示 积分 希 处 于 上 饱和 区 ， 
> 输出 为 0， 表示 积分 顺 处 于 正 背 范围 之 内 ， 
> 输出 为 一 1， 表 示 积 分 器 处 于 下 饱和 区 。 

















Clock 
Constant Sum 
20 
和 
说 明 ， 有 
(1) 积分 器 的 初始 值 为 0; U 
人 站 
(2) 【External reset】 的 设置 
为 : falling。 LI 本 
0 1 2 3 4 5 2: 了 8 9 10 


如 果 将 初始 值 设 置 为 5， 则 结果 为 : 


25 


20 














35 


30 


(1) 积分 器 的 初始 值 为 5; 
(2) 【External reset】 的 设置 为 :rising 。 














25 


20 


(1) 积分 器 的 初始 值 为 5; 
(2) 【External reset】 的 设置 为 ， either。 














(1) 积分 器 的 初始 值 为 5; 
(2) 【External reset】 的 设置 为 ， level。 











连续 系统 仿真 《采用 积分 模块 ) 








弹性 恢复 力 “及 
C 阻尼 力 CX 
系统 动力 方程 11IX 十 CY 十 KX 三 0 
系统 参数 : 1101. 三 0 全 2 
Cs 
对 原动力 方程 作 变换 : X 三 一 一 X 一 一 X 
1 711 
代入 参数 : 2 


假定 有 如 下 初始 条 件 ，”X(0) =1 ， 0O) =0 





惯性 力 


初始 值 为 1 





4 加 
Ce velocity displacement Scope 
NS 和 


177 177 Gain1 





和 
X 三 一 0.2X 一 0.4x7 


ROUTERAOEO II 


和 











-0.8 
0 


连续 系统 仿真 《采用 传递 函数 模块 ) 


WE 
攻 1IX ”惯性 力 
六 弹性 恢复 力 避 人 
系统 动力 方程 : 1 
经 过 Laplace 变换 ， 有 :1118“X (8S)TcsX(C)+TKXK(C) = 开 (9) 
传递 函数 为 : 

GO - X(C) 1 1771 02 


FG) ms +cs+ 2 CC 大 +0.2s+0.4 














1 711 
已 知 条 件 : 假设 系统 处 于 零 平 衡 位 置 ， 外 力 函 数 为 幅 值 为 1 的 阶 跃 函数 。 : 






Step Gain2 


本 国 
11] 十 [ 光 人 E 万 和 velocity displacement 






EEC 





X =-0U.2X 一 0.4xX 二 0.27 人 
X (9) 1 采用 传递 函数 模块 所 建 的 模型 
Cr(y) = 三 TREE 
(8S) 118 十 C8 十 大 和 
1 /1717 本 
S +0.2s+0.4 
和 十 S 十 Step Transfer Fcn Scope 
1 7111 | 
0.2 
汪 分 子 输入 : 0.2 分 母 输 入 : [1，0.2，0.4] 


whowasl 示 

















se (hd 














采用 幅 值 为 1 的 阶 跃 函数 激励 ， 随 着 时 间 的 推移 ， 速 度 和 加 速度 项 逐渐 趋 于 0， 
从 而 达到 静 平 衡 位 置 ，-0.4x 十 0.2*1 一 0， 即 Xx=0.5 


癌 量 线性 系统 


信号 线 传输 的 信号 可 以 是 标量 形式 ， 也 可 以 是 向 量 信号 形式 。 


问 量 信和 号 线 


在 Signals & Systems 模块 库 中 有 一 个 
Mux 模块 ， 它 可 以 将 多 个 标量 信号 合成 为 一 个 


四 量 信号 ， 输 出 为 列 同 量 。 


与 Mux 相对 应 的 模块 是 Demux 模块 ， 它 
用 来 将 回 量 信号 分 解 成 一 组 标量 信和 号 。 


在 Mux 模块 和 Demux 模块 中 必须 正确 指定 输入 或 输出 的 个 数 。 下 图 是 
将 三 个 标量 信号 合成 为 一 个 回 量 信号 ， 然 后 再 分 解 为 三 个 标量 信和 号 。 


输出 为 : 国 在 Gain 中 输入 [1,， 1]， 在 Multiplication 中 


SCFEE 全 症 
四 


一 [本 4 4 4 


Displa 
Constant1 Gain 的 


若 在 增益 模块 中 选择 Matrix (K * 由， 则 结果 如 下 。 


] 
输出 为 : 四 在 Gain 中 输入 [1，1]， 在 Multiplication 中 选择 Matrix (K*u) 
B 1 
ae 一 II 
四 天 Display 才 


Constant1 Gain 


答 出 为， | 在 Gain 中 输入 2， 在 Multiplication 中 选择 


-下 Element-wise (K.*UuU) 

了 | 一 由， 

Constant 一 一 > 5] 1 X2 一 : 
Display 


Constant1 Gain 


状态 空间 模块 
单 目 由 度 系统 的 单位 脉冲 啊 应 


11IX 十 CXY 二 AKC=O(f) 
3 十 X 直 2 三 (及 


网 辐 - 全 和 国 了 
区 | + lmD=47+BaOW) 
元 il lg 光 ie 





答 出 量 : 民 三 CI 要 求 输出 位 移 量 ， 则 : C=[ 0| 
近似 单位 脉冲 为 两 个 有 微小 时 间 差 的 阶 路 信号 的 差 : 
5() =100x(D) -100x4 一 0.0]) 


一 四 个 选项 : 0.01，100，0,， 0 
天 


丰 Step 
C 十 X' = AX+BUu 
三 y = CX+Du 
SS 


State-Space Scope 














11IX 十 CX 二 AX=O(1) sepi 一 四 个 选项 ,0， 100, 0,， 0 
NS | 色 
0.3 
有 和 全 全 环 
三 二 机 本 下 仆 和 生生 扑 生 相 和 1 1 L_ - -| 
0 人 
ee 人 下 - 
和 - 
Ze CK C=[L 0|] 
1 RE 5 - 
5(D) =100x(D)-100x0 -0.01) ee ] 
四 -| 
本 和 - 
四 0 了 
02 
本 10 15 20 25 30 35 40 


非 线性 系统 的 模拟 


例子 1: 小 于 由 两 个 喷射 式 发 动机 推动 在 光 请 的 平面 内 运动 。 
各 小 车 的 速度 和 位 移 之 和 为 负 值 ， 则 局 动 堪 边 的 发 动机 ， 奉 小 
车 的 速度 和 位 移 之 和 为 正 值 ， 则 局 动 右边 的 发 动机 。 欣 制 的 目 
标 使 小 车 静止 在 原点 。 此 模型 类 似 于 卫星 的 位 置 控制 过 程 。 





必 动力 学 方程 : 17IX 王八 信 苇 下 ss 
171 


ss X 二 X<I(O 
层 设 ， 有 亚 = 一 |]，X+X>0O 
已 = 一 10 


oa 付 





111 = 9Kg ， 


x 动力 学 方程 ， ME= 已 宗 =0.2F 


0 : 


| 











一 X 二 X<IO 
ax 
三 才 


ES 


假设 ， 1 =3Kkg， 上 =1 一 10 


系统 模型 初始 位 移 为 十 1 


假定 小 车 初始 项 
止 ， 并 且 位 移 量 
下 为 十 1 
Sign 1/m Velocity| Displacement 


符号 函数 模块 说 明 


ee | -XY 一 X>0 加 四 
| ee 时 aa 遇 Display 
L Sign ”Display | -5| 
日 Constanti 


| 


Constant1 


系统 模型 初始 位 移 为 十 1 










假定 小 车 初始 静 


Sign 1/m Velocity| Displacement 为 可 司 1 


用 一 个 二 维 图 形 模块 来 绘制 仿真 过 程 的 相 图 。 相 图 是 速度 相对 于 位 
移 的 变化 图 。 为 了 观测 时 间 ， 给 模型 加 一 个 时 钟 模块 。 





暗 关 辆 帮 司 


Sign 1/m Velocity| Displacement 


Display 











注 : 在 MATLAB 中 使 用 plot (xout(:,1)，xout(:,2)) 画 出 的 图 形 ， 和 XY Graph 的 结果 一 致 。 


为 了 得 知 小 车 何 时 到 达 原 扣 并 使 仿真 停止 ， 对 模型 增加 一 个 
网 辑 判断 ， 以 使 目标 达到 时 仿真 目 动 停止。 

对 于 本 算 例 ， 可 以 认为 小 车 的 位 移 和 速度 的 绝对 值 之 和 小 于 
一 个 较 小 量 ， 如 0.01 时 ， 观 认为 目标 已 经 达到 并 结束 仿真 : 


思 十 四 < 0U.01 
二 证 
Sign 1/m Velocity Displacement 
xYGraph 
Clock Display 


g | 


Relational Stop Simulation 
Operator 


Abs1 Constant 


如 果 想 得 到 小 车 位 移 的 响应 时 程 ， 可 在 模型 中 加 入 示波器 
模块 ， 如 下 图 所 示 : 












Velocity Displacement 
7. 图 
xYGraph 


Display 


-着 


Scope 


Clock 


Relational Stop Simulation 
Operator 


Abs1 Constant 


仿真 结果 如 下 页 图 形 所 示 。 










Sign 1/m 


Velocity| Displacernent Scope 0.8 


性 














-=-=--r--= TsTTTTTT 1 人 
| | 
Gek 和 
四 下 二 E 上 
Abs | 2 将 
CarSteStop Simulation ee | 
Abs1 Constant 1 1 2 5 4 。 6 8 





该 图 形 是 采用 MATLAB 环境 下 的 plotltout, xout(:,1)) 男 出 的 ， 
结果 与 示 波 占 的 结果 一 致 。 由 图 形 可 看 出 ， 仿 真 在 7.6 s 时 停 
止 ， 这 与 模型 框图 中 的 时 钟 显 示 模 块 Display 的 结果 一 致 。 小 车 
系统 在 7.6 s 时 达到 平衡 。 


中 极 跳 系统 的 动态 仿真 


蹦极 跳 时 一 种 挑战 号 体 极限 的 运动 ， 蹦 极 者 系 看 一 根 弹性 绳 从 忆 
处 的 桥 染 《或 山崖 等 ) 癌 下 跳 。 在 下 落 的 过 程 中 ， 咒 极 者 几乎 处 于 失 
重 状态 。 按 照 牛 顿 运动 规律 ， 目 由 下 落 的 物体 由 下 式 确 定 : 


11IX 三 1118 一 QIX 一 0 ee 


mm 为 人 体 的 质量 9 为 重力 加 速度 
位 置 x 的 基准 为 桥 染 的 基准 面 


如 果 人 体系 在 一 个 弹性 常数 为 K 的 弹性 绳索 
上 ， 定 义 绳索 下 中 的 初始 位 置 为 0， 则 其 对 落 
体位 置 的 影响 为 : 


一 KK，X>0 
D(X) = 三 
(U， X<(O 





因此 整个 蹦极 系统 的 数学 模型 为 : 
KX， X>O 


111X 一 1118 +D(z) 一 G 无 -Ga D(X) 有 1 


典型 的 具有 连续 状态 的 非 线性 系统 
设 桥架 距离 地 面 为 50 m， 即 m2=50 

蹦极 者 起 始 速度 为 0 
蹦极 者 的 起 始 位 置 为 绳索 的 长 度 30 m， 即 由 =30 0 





其 余 的 参数 : mw 
K 一 20， a 王 ai 一 1;， mm 一 70 kg，g9 王 10 m/s: 





初始 条 件 :  xX(0) =-30; X(O)=0 


1 人 之 二 118 十 已 (2 一 G 天 一 Go 让 x(O)=-30，F#(OD) =0 
关 = 7 [mg +DOD) 一 @ 无 一 ao] 厂区，x>0 
0， 人 


初始 条 件 0 初始 条 件 -30 
全 便 









Fcn 二 
spring 
constant 


Switch 人 This output the distance 


to the ground of jumper 


说 明 : 〈1) 设置 仿真 时 间 0 一 100s。 为 了 使 曲线 光滑 ， 可 设置 最 大 仿真 步 长 为 0.1。 
(2) 在 MATLAB 环境 下 输入 m=70; g=10; k=20; a1=1; a2=1; 


双击 输入 


Switch 模块 介绍 : Threshold 
Input 1 ee 中 
Input 2 四、 output 
Input 3 
Switch 


Pass through input 1 when input 2 is greater than Or equal to 
threshold; otherwise, pass through input 3. The inputs are numbered 


top to bottom (or left to right) 








国 Threshold 默认 为 0 攻 
Constant Constant 

国 码 广 一 | 回 
Constant2 Switch Display Constant2 Switch Display 


Constant1 


Constant1 


























| | | | | | | | | | | | | | | 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
有 和 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
FE==E==Eesk--u=ddccd- 站 下 村 人 生生 和 全 入 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 区 1 1 1 1 1 1 1 1 1 
区 F--r--r-tr--r--r--+--+--+--+--JR | 
\ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 了 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 攻 和 1 1 1 1 1 1 1 1 1 
上 和 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
A 户 1 1 3 1 1 1 1 1 1 1 区 1 1 1 1 1 1 1 1 1 
| 1 1 1 1 1 1 1 1 1 去 \ 1 1 1 1 1 1 1 1 1 
入 本 本 
1 1 1 1 1 1 1 1 正 1 1 1 1 1 1 1 1 
直 1 1 1 1 1 1 1 1 1 二 1 1 1 1 1 1 1 1 1 
一 1 1 1 1 1 1 1 1 1 后 1 1 1 1 1 1 1 1 1 
Ti En 
】 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 中 1 1 1 1 1 1 1 1 1 
名. 4- 医 | 
T T T T T 吕 [@ 罩 | | | T T T T 
[四 |】 1 1 1 T 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 四】 1 1 1 1 1 1 1 1 1 
操 人 本 
| 1 | T T T T T SU 扣 1 1 | | i T T f T 
(1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T 1 1 
1 1 1 1 1 1 1 1 1 CO 1 1 1 1 1 1 1 1 1 
| LT 1 1 1 1 1 | 
T 1 1 1 1 1 1 1 1 王 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 i T 1 1 1 1 1 1 1 1 1 1 T i 1 1 
上 上 上 上 | 上 | 上 口 | | 上 | | | | 上 
避 吕 总 口 所 避 局 口 局 避 吕 = 全 吕 吕 口 吕 局 口 它 本 局 
区 吕 豆 ) GU 一 二 G 中 吕 区 加 re) 于 G GU 去 二 忆 


Scope 
站 
十 
This output the distance 
to the ground of Jumper 









会 


spring 
constant 


初始 条 件 0 ”初始 条 件 -30 






Air resistance 





吉 在 MATLAB 环境 下 键入 : whos 


可 得 到 : 
whos 
Name Size Bytes Class 
al 1X1 8 double array 
a2 1X1 8 double array 
g 1X1 8 double array 
k 1X1 8 double array 
m 1X1 8 double array 
tout 1000X1 8000 double array 
Xxout 1000Xx2 16000 double array 


tout 为 时 间 序 列 ，xout 第 一 列 为 位 移 时 间 序 列 ， 第 二 列 为 速度 时 间 序 列 。 阁 在 
MATLAB 环境 下 键入 : plot ( tout，xout ( :，1) )， 将 得 到 示波器 Scope 显示 的 结果 


如 果 想 在 MATLAB 环境 下 使 用 plot 命令 画 出 示波器 Scope1 显 
示 的 图 形 ， 可 在 系统 模型 中 加 入 out 模块 。 


mg 


Constant 








Position Scope 


onstanti 
spring 
a2*abs(U)* constant 
Fcn 
Out1 


This output the distance 





Velocity 





Air resistance 


[ LO | 


to the ground of jumper 


此 时 在 MATLAB 环境 下 键入 : whos 


可 得 到 : 
whos 
Name Size Bytes Class 
al 1X1 8 double array 
a2 1X1 8 double array 
g 1X1 8 double array 
k 1X1 8 double array 
m 1X1 8 double array 
tout 1000x1 8000 double array 
Xxout 1000X2 16000 double array 
yout 1000xX1 8000 double array 


其 中 yout 为 out 模块 的 输出 结果 。 此 时 若 在 MATLAB 环境 下 键入 : 
plot ( tout，yout ) 


将 得 到 系统 模型 中 示波器 Scope1 显示 的 结果 。 





constant 





This output the distance 
to the ground of jumper 


对 于 体重 为 70 kg 的 蹦 
极 者 ， 此 系统 是 不 安全 
的 ， 因 为 咒 极 者 与 地 面 
之 间 的 距离 出 现 了 负 值 





Scope 显示 的 结果 

















Scope1 显示 的 结果 











SIMULINK (0) 


子 系统 及 其 封装 技术 


本 章 学 习 内 容 和 目的 


@ 掌握 Simulink 子 系统 的 建 模 方法 
@ 掌握 Simulink 子 系统 的 封装 技术 


在 前 面 的 章节 中 ， 介 绍 了 使 用 Simulink 进行 建 模 的 基 
本 方法 。 使 用 这 些 方法 基本 可 以 创建 任何 物理 系统 的 模型 。 
然而 随 看 系统 越 来 越 复 杂 ， 用 这 些 基 本 操作 创建 的 
Simulink 模型 变 得 越 来 越 庞大 而 难于 读 介 。 在 本 章 中 ， 将 
介绍 一 系列 的 Simulink 的 特殊 处 理 技术 ， 使 得 模型 变 得 更 
加 简捷 和 易 懂 易 用 。 


本 章 首 先 介绍 一 种 类 似 于 程序 设计 语言 中 的 子 程序 的 
处 理 方法 一 Simulink 子 系统 ， 然 后 讲解 一 种 更 加 好 用 的 
封 厂 子 系统 的 技术 。 


Simulink 子 系统 


绝 大 多 数 程序 设计 语言 都 有 使 用 子 程序 的 功能 ， 例 如 MATLAB 的 子 
程序 一 一 M 文 件 。 随 独 模 型 变 得 越 来 越 大 、 越 来 越 复 杂 ， 人 们 很 难 轻易 读 
介 它 们 。 在 这 种 情况 下 ， 子 系统 通过 把 大 的 模型 分 割 成 几 个 小 的 模型 以 使 
得 整个 模型 更 加 简捷 、 可 读 性 更 高 ， 而 且 这 种 操作 并 不 复杂 。 


创建 Simulink 子 系统 有 两 种 方法 : 


(1) 对 已 经 存在 的 模型 的 某 些 部 分 或 全 部 使 用 沈 单 命 令 【Edit>Creat 
Subsystem 将 你 性 压缩 转换 ， 使 之 成 为 子 系统 ; 


(2) 使 用 Subsystems 模块 库 中 的 Subsystem 模块 直接 创建 子 系统 。 


Simulink 子 系统 的 两 种 作用 : 
(1) 系统 模型 更 加 简捷 和 可 读 性 高 ; 
(2) 子 系统 可 以 反复 调用 ， 节 省 建 模 时 间 。 


压缩 子 系统 


以 一 个 例子 说 明 压 缩 子 系统 的 使 用 方法 。 例 如 : 
X 十 U.4X 二 0.8X 一 SIn(T) 


方程 可 转化 为 : 
X 三 -0.4cx 一 0.8X 二 SIn(T) 


是 | Displacemeht “ Scope 


芯 


Sine Wave 





芯 


操作 步骤 : 


(1) 使 用 范围 框 将 要 压缩 的 子 系统 的 部 分 选中 ， 包 括 木 块 和 信 
号 线 ;， “注意 : 只 能 使 用 范围 框 ， 而 不 能 使 用 Shift 逐 个 
选 定 ) 

(2) 在 模块 窗口 选项 中 选择 【Edit>Creat Subsystem 】， 
Simulink 将 会 用 一 个 子 系统 模块 代 蔡 被 选中 的 模块 组 ; 


《3) 进行 模型 美观 调整 。 


Eee 口 和 输出 端口 


四 王 Wave 
上 二 Wave 国 
加 


Subsystem Subsystem 


压缩 子 系统 后 的 系统 模型 美观 处 理 后 的 系统 模型 


Sine Wave Scope 


Subsystem 


| 双击 子 系统 ， 则 会 出 现 一 个 显示 
子 系统 内 容 的 新 窗口 。 






Displacement Quti 


Gain1i 


在 新 窗口 中 ， 除 了 原始 的 模块 外 ，Simuink 自动 添加 了 输入 模块 和 输出 模块 ， 分 
别 代表 子 系统 的 输入 端口 和 输出 端口 。 


两 点 说 明 : 

(1) 子 系统 窗口 无 需 保 存 ， 只 需 保 存 主 程序 出 口 即 可 。 保 存 
主 程序 窗口 后 ， 子 系统 窗口 目 动 得 以 保存 ; 

(2) 菜单 命令 【Edit>Creat Subsystem】 没有 相反 的 操 

作 命令 ， 也 就 是 说 ， 一 旦 一 组 模块 压缩 成 了 子 系统 ， 就 没有 

可 以 直接 还 原 的 处 理 方 法 (undo 除外 ) 。 因 此 ， 一 个 理想 


和 得 着 和 生 人 先 将 模型 进行 保存 ， 作 为 
份 。 


以 上 例子 的 运行 结果 
假定 系统 有 如 下 初始 条 件 : X(O) = XGO)=(0 


要 求 采 用 原始 模型 “没有 压缩 子 系统 的 模型 ) 和 压缩 子 系统 
的 模型 进行 仿真 ， 仿 真 时 间 30s。 对 比 结 宋 。 


压缩 子 系统 的 模型 


| 司 。 同 e 







Sine Wave 


Sine Wave Scope 


Subsystem 








子 系统 模块 


绍 的 压缩 子 系统 的 方法 外 ， 也 可 以 直接 使 用 子 系统 模块 ， 在 子 
系统 窗口 中 进行 创建 。 


要 使 用 子 系统 模块 创建 子 系统 ， 先 从 Subsystems 模块 
库 中 拖 电 一 个 子 系统 模 块 到 模型 窗口 中 ， 然 后 双击 子 系统 模 
块 ， 就 会 出 现 一 个 子 系统 的 编辑 窗口 ， 子 系统 的 建立 可 以 在 该 
窗口 中 进行 建立 。 


例子 : 模拟 下 图 所 示 的 弹 赞 一 质量 系统 的 运动 状态 。 





2 > 
ZX ZX TD 人 TAR 


本 1 着 十 大 5 十 所 (一 六 )=0 
由 > > 
一 m -一 | 
| 
KR 171 


112 思 一 局 一 六 ) 二 后 一 为) = 


.. ] 
二 nl 


10 入 一 局 Op 一 为)=0 
[Co 一) 





| 
En 
/1 


.1 
单个 质量 的 运动 方程 : -oo 本 [ 0 一 志 ) | CC 一 ) 


双击 后 可 显示 出 第 一 个 质 










”| 量 的 子 系统 图 ， 下 图 所 示 
于 由 
第 二 和 第 三 个 质量 子 系统 中 ， 
三 个 增 义 模块 分 别 输入 : 
Right Spring-mass3 一 
free 





第 三 个 质量 : 上 K。， 0， ms 


人 
Left 
-证 spring 


( ) | 注意 ;在 子 系统 的 信 
pisp “| 号 输入 端 要 使 用 一 个 





国 
十 Resultant 


Right ”Sum1 必 spring 


区 时 Positio 


输入 模块 ， 在 信号 输 
出 器 要 使 用 一 个 输出 
模块 。 






Block displacement 


假定 系统 的 初始 条 件 为 : 第 一 个 质量 上 有 初始 位 移 xi(0)=1， 则 在 第 
一 个 质量 子 系统 框图 中 的 位 移 模 块 中 输入 初始 条 件 1。 在 MATLAB 环境 





下 输入 : 
K1=1; KK2=2; kk3=4; m1=1; m2=3; m3=2 

系统 仿真 结果 如 下 : 
生 说 明 : 若 想 输 出 其 它 两 
人 DR 个 质量 的 轨迹 ， 可 在 模 
ET 3 型 中 增添 示波器 模块 
四 -上 Scope ， 或 增添 输出 模 
FE 块 Out 即 可 。 











封 芍 模块 


封装 技术 是 将 Simulink 子 系统 “包装 ”成 一 个 模块 ， 并 可 
以 如 同 使 用 Simulink 内 部 模块 一 样 使 用 的 一 种 技术 。 每 个 封 
闭 模 块 都 可 以 有 一 个 自 定 义 的 图 标 和 一 个 用 来 设 定 参 数 的 对 话 
框 。 参 数 的 设 定 方法 也 与 Simulink 模块 库 中 的 内 部 模块 完全 
相同 。 


作业 : 


封 闻 踊 极 系统 。 要 求 : 封装 后 的 蹦极 子 系统 只 有 一 个 
和 输出 端口 ， 封 厂 后 子 系统 的 参数 设置 包括 中 极 者 的 体重 、 
弹性 绳索 的 弹性 单数 。 通 过 仿真 分 析 蹦 极 系统 在 下 述 情况 
下 是 否 安全 ， 并 绘制 啊 应 的 啊 应 曲线 ; 


(1) 蹦极 者 体重 80 kg， 弹 性 绳索 的 弹性 常数 为 30; 
(2) 蹦极 者 体重 70 kg， 弹 性 绳 款 的 弹性 利 数 为 20。 


